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Rapid development in electrical technology has imposed strong challenges 
to modern power system. Power quality has become a great concern due to 
proliferation of power electronic technology in modern electrical loads. 
Specifically for lighting load such as compact fluorescent lamps (CFLs), one 
of the concerning issues is harmonics. CFL is a cost-competitive and energy 
efficient compared to incandescent lamp. Inevitably, CFL produces 
harmonics current due to nonlinearity behaviour of the electronic ballast 
circuit. This paper presents a study on the widespread installation of CFL 
lamps in electrical power network. Initially, the harmonic current 
characteristics of local-branded CFL was identified from laboratory 
measurement. Then, a simulated CFL model was developed in 
MATLAB/Simulink to replicate the identified characteristics. The same step 
was repeated for other two different brands where eventually all models 
were embedded into a distribution network. The results show that at low 


voltage level, with installation more than 50 units for each type of CFL, the 
harmonic voltage distortion exceeded the 8% total harmonic distortion 
(THD) limit as stipulated in EN50160 standard. However, at higher voltage, 
the amount of THD decreased to average 0.94% and further down to average 
0.28% at small transmission voltage level. 
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1. INTRODUCTION 

The adoption of power electronic based technology load is widely spreading and continuously 
increasing. The technology applications vary from as small as residential lightings to heavy industrial drives. 
For many years, fluorescent lamp (FL) has received wide acceptance in residential lighting replacing 
conventional incandescent lamp because of higher efficiency and longer lifespan [1]. With advancement in 
power electronic technology, the fluorescent lamp design has also evolved and leads to the invention of 
compact fluorescent lamp (CFL). CFLs are still widely available in the market because of efficiency and 
economic value [2]. The global lighting sales report released by international energy agency (IEA) shows 
CFL holds 28% of market share in 2020 [3]. 

The electronic ballast of CFL has nonlinear behaviour, hence it distorts the sinusoidal signal of AC 
current [4]. Distorted sinusoidal signal contains harmonics as shown in Figure 1. In the figure, the reference 
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signal is the fundamental sine wave at 50 Hz frequency while the harmonics are others that exist at multiples 


order of fundamental frequency which are 150 Hz (3 order), 250 Hz (5" order), and 350 Hz (7" order) sine 
waves accordingly. 


50Hz fundamental 


L e freq. signal 
yo Distorted signal 


5th harmonic 


signal ` 3rd harmonic signal 


7th 
` harmonic ¢ 
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Figure 1. Decomposition of distorted sinusoidal signal 


When the distorted current wave passes through impedance in the electrical network it further 
induces distortion in the voltage waveform that causes adverse power quality issues [5], [6]. Among serious 
problems caused by the harmonic currents in the power system are power loss, overheating of conductors, 
overloading of motors and transformers; undesired tripping of protection devices, malfunction of logic 
controllers and poor power factor [7]-[9]. The operation of nonlinear loads such like switch mode power 
supplies (SMPS), power converters, variable speed drives and battery chargers continues to increase within 
the power system. Research by Kalair et al. [10] claims reduction in power factor around 60% and an 
increment in line losses more than 2% due to widespread use of nonlinear loads. Thus, with proliferation in 
nonlinear loads and widespread installation of energy saving alternatives like CFLs and LEDs, it is important 
to model the complex interactions and to analyze the impacts on electrical network. 

This paper focuses on modelling harmonics current injection of CFL using MATLAB/Simulink to 
evaluate the effect of CFL widespread installation on power network. Three different units of CFL are 
studied to include possible varieties in harmonic characteristics. The actual data of current measurements are 
used as reference to model the harmonic current sources. The interaction of multiple CFLs at different scales 
of installation is simulated using the proposed harmonic current model with a theoretical electrical network. 
The simulation model is used to assess the total harmonic distortion (THD) indices in the network at various 
observation nodes with relate to size of installation. 

The remaining part of the paper is as follows. Section 2 discusses some scientific literature related to 
nonlinear load models to characterize the harmonics current injection. The MATLAB/Simulink model for 
CFL and descriptions of the study case is given in section 3. Section 4 shows the results and analysis. Lastly, 
section 5 concludes this study. 


2. RELATED WORK 

In general, circuit-based and measurement-based are two main methods used in harmonic 
modelling. Analyzing harmonic models with simulation softwares like power system computer aided design 
(PSCAD)/electro magnetic transient DC (EMTDC) and MATLAB offer fast and reliable solution. Circuit- 
based models have been used to assess the harmonic injections from nonlinear loads in several researcher 
works. Such as in [11]-[14] modelling the equivalent circuit of full wave rectifiers in MATLAB/Simulink to 
characterize harmonic components. In [15], [16] using a simplified model of SMPS to study harmonics 
propagation and interaction among end-user loads in PSCAD and [17] using estimated equivalent parameter 
values of SMPS in modern household appliances for aggregated harmonic analysis. The diversity of 
harmonics spectrum generated with different set of parameter values in the CFL equivalent circuit was 
studied by [18] using alternative transients program-electromagnetic transients program (ATP-EMTP) and 
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MATLAB programs. The study suggested that the behavior of distorted current is not predictable, thus, is 
vital to find the parameters that describe the distorted current. 

For measurement-based model studies the nonlinear load has always been modelled as a source of 
voltage or current depending on the purpose of study [19], [20]. Research by Skamyin et al. [21] deals with 
deploying passive filter in electrical network in the presence of harmonic currents. The equivalent nonlinear 
load is represented in the form of current sources at respected harmonic frequencies. The summation or 
subtraction of the currents due to passive filter deployment is further used in determining the efficiency of the 
filter. According to Watson et al. [22] the spectra of harmonics current from different CFLs were determined 
by laboratory measurements resulting with representation of Norton equivalent circuit of all CFLs. Sola and 
Salichs [23] modeled harmonics generated by different CFL types with data from harmonic current 
measurements. The internal electronics circuits are considered as a “black box”. The harmonics spectra were 
modeled using Pspice simulation, analytically processed and eventually characterized by means of parametric 
conductance. A similar work observed by Sainz et al. [24] embraces the “black-box” concept from [23] with 
more comprehensive way of characterizing the waveform using nonlinear least squares method for estimation 
of parametric conductance. 

In real practice, the challenge of circuit-based method is to correctly attain the complete circuitry or 
circuit parameters, which neither is easy to obtain in modern appliances. On the other hand, several reliable 
measurement modelling methods have also been adopted. The differences are in the way measurements data 
being represented in harmonic model which could lead to a compromise between complexity and accuracy. 
Finally, in modelling any nonlinear load, rather than the modelling approach itself is important, the 
convenience of applicability in commercial simulation programs is also worth to be considered. 


3. RESEARCH METHOD AND MATERIAL 
3.1. Experiment setup 

The proposed study is accomplished through experimental and simulation work. The equipment 
setup for harmonics measurement in laboratory is shown in Figure 2. The CFL is connected to a circuit that 
feeds 240 Vrms and 50 Hz supply. The criteria of tested CFL units include identical brand but different color 
temperature, different brand but identical power rating and identical brand but different power rating. The 
rated values of the CFL lamps used in the experiment are given in Table 1. Fluke power quality analyzer 
(PQA) model 437-II and fluke AC current clamp i400 are used in measuring the harmonic currents. The fluke 
PQA features 16-bit resolution and 200 kS/s sampling rate whereas the current clamp sensitivity is 1 mA/A. 
Comparatively the sampling rate of the utilized PQA is a bit low which sets some limitation on the 
measurement results. All collected data is then processed using fluke power log (version 5.9) software 
installed in personal computer. 


Supply voltage Compact fluorescent lamp 


240 Vrms @ 50 Hz 


@& F Í 5 
Current 
probe 


Power quality analyzer Personal computer 


Figure 2. Experiment setup for CFL measurement 


Table 1. Rated values for experimental lamps 


Lamp brand Power (W) Voltage (V) Frequency (Hz) Colour 
Type 1 (brand A) 18 220-240 50/60 Cool day light 
Type 2 (brand B) 18 220-240 50/60 Cool day light 
Type 3 (brand B) 20 220-240 50/60 Warm white 
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3.2. Simulink model 

In this work, the proposed harmonic simulation model is developed using MATLAB R2021a with 
Intel 15-4590 CPU and 8 GB RAM computer. Using common Simulink block, two simulation models were 
developed to represent the electrical system, namely the CFL current model and electrical distribution model. 
For CFL current, there are three sub-models representing three different types of CFL lamp and for 
distribution model, a small theoretical network is considered. 


3.2.1. CFL harmonic current model 

The behaviour of the measured CFL harmonic current is replicated in MATLAB/Simulink using 
Fourier series function. The function simply transforms the fundamental and harmonics of current 
measurements in time domain into frequency domain. Such approach commences from the concept that any 
function waveform, S(t) can be transformed into fourier series in amplitude-phase form as given by [25]. 


S(t) = + DNL, [Ancos (=) - on] (1) 
where 7 is the the mean value of function, A, is the n™ harmonic amplitude, y, is the n™ harmonic phase 
shift, n is an integer index which is also the number of cycles the n™ harmonic makes in the function period, 
P is the function period, t is the time index and N is the number of observation terms. Thus, by identifying 
individually the magnitude and frequency of the measured data, the harmonics waveform is reconstructed by 
combining a sum of sinusoidal function at identified harmonic frequencies. The CFL model is represented 
using the Simulink function block and sinusoidal block with each phase has different functional block 
according to the type of CFL as shown in Figure 3. 
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Figure 3. CFL harmonics current model 


3.2.2. Theoretical electrical network model 

The electrical network is a large system comprises of numerous network components and 
configuration. For accurate representation of a specific network, adequate information of the components is 
necessary. Since the purpose of the modelling is to comprehend the idea behind the interaction of CFL 
harmonics with network impedance, therefore any typical small network is satisfactorily applicable. To 
purely assess the harmonic effects from CFL, the load system is limited only to a mix of three CFLs without 
other nonlinear loads. The CFL loads are modeled as current sources that inject various harmonic spectrums 
into the network. The severity of harmonic impacts is assessed in terms of THD in current (THD;) and THD 
in voltage (THD,,). The expressions of THD; and THD, are given in (2) and (3) respectively [26]: 


[Ee pital? 
THD;(%) = = x 100 


Il 


(2) 


where h is the order of harmonic, J, is the rms value of current component at frequency hf, and J, is the rms 
value of fundamental component of the current. Equivalent expression for THD, is: 


TAE 
THD,(%) = x 100 


ue (3) 
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where h is the order of harmonic, V, is the rms value of voltage component at frequency hf, and V; is the rms 
value of fundamental component of the voltage. The steps taken for investigating the CFL-network 
interaction are summarized as follows: 
a. A CFL-clustered load comprises of three different CFLs connected to low voltage network. Initial load of 
10 clusters is chosen and THD; is recorded at PCC4 which is the nearest point to the load. 
b. Increase load in a step size of 10 clusters per increment and repeat the step until total load reaches 50 
clusters. 
c. The THD, of the network is also necessary to be observed. The load size increment procedure is repeated 
and THD, is recorded at all PCCs for each event. 
The proposed model comprises of a three-phase supply network with illustrated size of CFL loads 
shown in Figure 4. 
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PCC 3 Transformer 
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Figure 4. CFL harmonic model in electrical network 


4. RESULTS AND DISCUSSION 

The current waveforms for CFL measured from the experiment in are presented in Figure 5. All 
waveforms almost similar in shape with small differences in magnitude. Some oscillations in the waveforms 
are due to measurement noise. The time-domain harmonic signals were analyzed using power log tool by 
fluke to be transformed into frequency-domain spectrum as shown in Figure 6. It shows that CFL type 3 
produces slightly higher current distortion under the same operating voltage. The spectrum considers 
harmonic order from 2™ to 20" because the effects of harmonics is usually prominent up to 25" order for 
lighting device; for harmonic order 21* and onwards the magnitudes are considerably small. 
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Figure 5. Current measurement waveforms for three different CFL types 
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Figure 6. Harmonic spectrum of current measurement for three different CFL types 


Based on the measurement results, the harmonic waveforms were reconstructed as shown in 
Figure 7. The waveforms in Figures 7(a), 7(b) and 7(c) represent harmonic current for CFL type 1, CFL type 
2 and CFL type 3 respectively. The shape of the modelled harmonic current approximately matches with 
measurement current. Some variations were probably due to measurement noise and the range of harmonics 
order. More important is the current amplitudes at lower frequencies are much higher particularly in third 
order and few other odd orders that give the sharp increase in current amplitude as observed in the waveform. 
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Figure 7. Harmonic waveforms of modelled current (a) CFL type 1, (b) CFL type 2, and (c) CFL type 3 


Table 2 summarizes the simulation results after increasing the load size. In general, the simulation 
results indicate the magnitude of harmonic current varies insignificantly with addition of harmonic sources, 
correspondingly the THD; remains constant around 78%. The THD; obtained is in line with experiment 
results achieved by [26] using three different CFL units whereby the THD; equals to 84%. It can be seen that 
increasing the cumulative harmonic disturbance by identical group of harmonic sources at the same 
proportion does not change the total harmonic current distortion at the same point of connection. 

On contrary, the total voltage distortion or THD, listed in the last column of Table 2 shows some 
increase with addition of harmonic sources. It can be seen in every cluster size, the THD, almost equals to 
arithmetic sum between the present and initial THD, values. This situation is consistent with the Ohm’s law: 


V=IxZ (4) 


where V is rms value of voltage, I is rms value of harmonic current and Z is impedance of network. 
Assuming Z is constant, the fundamental current and each harmonic current that pass through the impedance 
will cause a voltage drop across the impedance that eventually distort the waveform of voltage. Thus, it can 
be seen that the amount of voltage distortion at a location depends upon the size of the harmonic sources. 
Figure 8 illustrates the relation between various harmonic components and size of CFL loads. 


Table 2. Harmonic current magnitudes with varying size of CFL clusters 
No. of clusters %3" %5 %7” %9 %11 %13 %15 %THD; %THD, 


10 54.44 27.97 23.96 17.96 8.48 9.98 8.97 78.21 1.66 
20 54.43 27.96 23.96 17.96 8.48 9.97 8.96 78.19 3.31 
30 54.43 27.99 23.98 17.96 8.47 9.97 8.96 78.19 4.96 
40 54.43 27.96 23.95 17.95 8.47 9.97 8.96 78.18 6.62 
50 54.43 27.96 23.96 17.96 8.48 9.97 8.96 78.19 8.29 
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Figure 8. Relation between harmonic components and size of CFL loads 


Additionally, the interaction of harmonic currents with system impedance is further investigated. 
The observation of total voltage distortion, THD, is extended to other PCCs in the network. It is observed 
that the THD, at all PCCs shows increment trend with the increase in number of CFL loads. Figure 9 
illustrates the relation between harmonics voltage distortion and distance of common coupling point with 
different size of CFL loads. 
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Figure 9. Relation between THD, and distance of PCC 


Table 3 shows variation of harmonics distortion indices with respect to size of CFL loads at 
different coupling point. The obtained results show THD, in the network ranges from 0.06% to 8.29%, with 
the lowest distortion recorded at PCC1 and highest at PCC 4 for every cluster size. Evidently, the simulation 
results are justified with theoretical explanations based on Ohm’s law given in (4) as the greater impedance 
of the feeder, the higher the amount of voltage distortion. Since PCC 4 is the farthest point away from the 
main supply, so it accounts for the largest impedance in the network. Therefore the THD, recorded at this 
point with the maximum size of CFL cluster is the highest. Overall, it can be concluded that the amount of 
THD, is a function of load current, harmonic distorted currents and system impedance. The point of 
connection of harmonic-infested load in the network will determine the amount of distorted harmonic voltage 
at the chosen PCC. 

In this part, the THD; results were also investigated at PCC 1, PCC 2 and PCC 3 during the 
simulation. Clearly, the current harmonic distortion shows the highest magnitude at low voltage feeder which 
is PCC 4 where the harmonic currents initiated. On contrary, the lowest magnitude of current distortion is 
recorded at PCC 1. Primarily, this is due to the attenuation since PCC 1 is the most distant point from where 
the harmonic current sources are located. Using the same concept, it can easily clarify the reason for THD; in 
PCC 3 is higher than PCC 2. 

On additional note, Table 3 also indicates amplification in harmonic current distortion with 
increasing size of CFL loads for every PCC except for PCC 4 which is against the prior finding. It can be 
seen that the increment trend is not similar and unpredictable. A possible explanation to this variation is the 
effect of X/R ratio which is the ratio of the system reactance to resistance. In the network model, each PCC is 
identified with different level of voltage and accordingly its X/R ratio also has different value. The 
differences in X/R ratio will possibly cause diversity effect which harmonic currents with different phase 
angles are formed resulting with variation in harmonic spectrum. Eventually, it can be concluded the severity 
of harmonic current distortion is not purely dependent on the size of harmonic sources but also on attenuation 
factor and diversity factor of the network. 
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Table 3. THD, and THD; at all PCCs with varying CFL clusters 


PCC 1 PCC 2 PCC 3 PCC 4 
No of clusters %THDv %THDi %THDv %THDi %THDv %THDi %THDv %THDi 
10 0.06 3.67 0.07 16.03 0.18 41.00 1.66 78.21 
20 0.07 6.95 0.11 24.80 0.34 43.80 3.31 78.19 
30 0.09 9.88 0.16 29.99 0.51 44.78 4.96 78.19 
40 0.10 12.50 0.20 33.34 0.68 45.28 6.62 78.18 
50 0.13 14.85 0.25 35.66 0.85 45.59 8.29 78.19 


5. CONCLUSION 

It is found that increasing the generation of harmonic currents from purely similar group of 
harmonic source causes varying reaction to the network. The simulation results indicate the magnitude of 
harmonic current does not vary significantly with an increase in harmonic sources as the THD; remains at 
78% for the chosen point of observation. However, voltage distortion or THD, shows significant increment 
between 1.66% and 8.29% when the size of harmonic sources is increased five times larger. The severity of 
voltage distortion at other coupling points is directly determined by its equivalent impedance. The results 
show THD; varies between 3.67% and 78.19% while THD, varies between 0.06% and 8.29%, with the 
lowest recorded at PCC far from harmonic sources and the highest at PCC near to harmonic sources. It is 
clear, the installation of CFLs in large quantities profoundly reduces the power quality in terms of harmonic 
distortion especially at low voltage network. 
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